function f=fxyz(p)
x=p(1);y=p(2);z=p(3);
f=x+y^2/x/4+z^2/y+2/z;
